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Abstract 



We consider a critical branching particle system in M , composed of individuals of 
a finite number of types i E {!,..., K}. Each individual of type i moves independently 
according to a symmetric Oj-stable motion. We assume that the particle lifetimes and 
offspring distributions are type-dependent. Under the usual independence assumptions 
in branching systems, we prove extinction theorems in the following cases: (1) all the 
particle lifetimes have finite mean, or (2) there is a type whose lifetime distribution has 
^ ' heavy tail, and the other lifetimes have finite mean. We get a more complex dynamics 

by assuming in case (2) that the most mobile particle type corresponds to a finite-mean 
i^ . lifetime: in this case, local extinction of the population is determined by an interaction 

^— ^ , of the parameters (offspring variability, mobility, longevity) of the long-living type 

P^ I and those of the most mobile type. The proofs are based on a precise analysis of the 

occupation times of a related Markov renewal process, which is of independent interest. 
Keywords: Critical branching particle system; Extinction; Markov renewal process. 
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5^ ; 1 Introduction 

In critical branching and migrating populations, mobility of individuals counteracts the ten- 
dency to asymptotic local extinction caused by the clumping effect of the branching. In fact, 
convergence to a non-trivial equilibrium may occur in a spatially distributed population 
whose members perform migration and reproduction, even if the branching is critical, pro- 
vided that the mobility of individuals is strong enough. This behavior has been investigated 
in several branching models, including branching random walks [7, 9], Markov branching 
systems (both with monotype [6] and multitype [3, 4, 8] branching), and age-dependent 
branching systems [15]. 

In [15] Vatutin and Wakolbinger investigated a monotype branching model in Euclidean 
space W^, in which each particle moves according to a symmetric a-stable motion, and at 
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the end of its lifetime it leaves at its death site a random number of offsprings, with critical 
offspring generating function /(s) = s + |(1 — sy~^'^, f3 G (0, 1]. It turned out that, if the 
initial population is Poisson with uniform intensity and the particle lifetime distribution has 
finite mean, such process suffers local extinction ii d < a/P, while for d > a/f3 the system 
is persistent, i.e. preserves its intensity in the large time limit. This result is consistent with 
the intuitive meaning of the population parameters: the exponent a > is the mobility 
parameter of individuals in the sense that a smaller a means a more mobile migration (i.e. 
more spreading out of particles) which is clearly in favor of persistence; /3 is the offspring 
variability parameter, meaning that a smaller /3 causes a stronger clustering effect in the 
population, which favors local extinction due to criticality of the branching. If the lifetime 
distribution has a power tail t~^ for some 7 G (0,1], then the critical dimension is a'-f/p. 
Again, it is intuitively clear that long lifetimes (i.e. small 7) enhance the spreading out 
of individuals. However, Vatutin and Wakolbinger discovered that, in contrast with the 
case of finite- mean lifetimes, if the lifetimes have a general distribution of the above sort, 
the "critical" dimension does not necessarily pertains to the local extinction regime: when 
d = a'y/f3 persistence of the population is not excluded. 

Our aim in the present paper is to get a better understanding about how population 
characteristics such as mobility, offspring variability, and longevity of individuals determine 
the asymptotic local extinction of branching populations. In order to attain this we deal with 
a multitype system, where the most mobile migration (corresponding to the smallest a) and 
the life-time distribution with the heaviest tail, may correspond to different particle types. 
More precisely, we consider a branching population living in M.'^, constituted of particles of 
different types i G K := {!,..., K}. Each particle of type i moves according to a symmetric 
Qjj-stable motion until the end of its random lifetime, which has a non- arithmetic distribution 
function F,. Then it branches according to a multitype offspring distribution with generating 
function /i(s), s G [0, 1]^, i G K. The descendants appear where the parent individual died, 
and evolve independently in the same manner. The movements, lifetimes and branchings 
of particles are assumed to be independent; the only dependency in the system is that the 
offsprings start where the parent particle died. In addition, we assume that the process 
starts off at time from a Poisson random population, with a prescribed intensity measure, 
and that all particles at time have age 0. Let M = {'mij)f:^ denote the mean matrix of 
the multitype branching law, that is 

where 1 = (1, 1, . . . , 1) G M^. We assume that f (s) = (/i(s), . . . , /i^(s)) ^ Ms, and that M 
is an ergodic stochastic matrix. This implies that the branching is critical, i.e. the largest 
eigenvalue of M is 1. 

For the system described above, here we investigate parameter configurations under which 
the population becomes locally extinct in the large time run. We deal first with the case 
when all particle lifetimes have finite mean and prove that the process suffers local extinction 
ii d < a//?, where the mobility parameter a = mini<j<xai is the same as in the Markovian 
case [8], and the offspring variability parameter /3 G (0, 1] is determined by 

X — (v, 1 — f(l — ux)) ~ x^^'^L{x) as X — )■ 0, 



where v denotes the (normahzed) left eigenvector of the matrix M corresponding to the 
eigenvalue 1, and L is slowly varying at in the sense that \\m.x^Q L{\x) / L{x) = 1 for every 
A > 0. In a way, this case is similar to the one with exponentially distributed lifetimes. 

Next we assume that exactly one particle type is long-living, i.e. its lifetime distribution 
has a power tail decay t^"', 7 G (0, 1], while the other lifetime types have distributions with 
tails decaying not slower than At~^ for some rj > 1, A > 0. We consider two scenarios. In the 
first one we assume that the most mobile particle type is, at the same time, long-living, and 
we prove that extinction holds when d < a'-f/p. Then we proceed with the most interesting 
case: the most mobile particle type corresponds to a finite-mean lifetime. In this scenario, it 
turns out that local extinction of the population is determined by a complex interaction of 
the parameters (offspring variability, mobility, longevity) of the long-living type and those 
of the most mobile type. Assuming without loss of generality that type 1 is the long-living 
type, we prove that the systems suffers local extinction provided that d < dj^, where 
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The positive number '-/ai can be considered as the "effective mobility" parameter of the 
long-living type. If 70:1 is very close to a (so that 7/0; and l/ai are approximately the 
same), then d+ is also close to «i7//3 and to a/ (3. Moreover, for fixed a,ai and 7, the 
parameter d+ considered as a function of f3, is decreasing, which is consistent with previous 
known results. 

The proofs of our results rely on a precise asymptotic analysis for the occupation times in 
the branching particle system. Some of our techniques combine parts of the approaches used 
in [3] and [15] adapted to our model, however the adaptation to our case is far from being 
straightforward. In Section 2 we provide a family tree analysis which allows us to compare 
the occupation times of the particle system with the occupation times of an auxiliary Markov 
renewal process. 

Then, in Section 3, we carry out the asymptotic analysis mentioned above by investigating 
the occupation times of all types in the auxiliary renewal process, as well as the asymptotic 
number of renewals in large time-intervals. This is the mathematical core of the paper, and 
we think it is interesting on its own right. To achieve this, we need to control the tail decay 
of the renewal times of all types simultaneously, which we were able to do assuming that 
there is only one long-living particle type. Therefore, in its present form our approach is not 
yet applicable to treat a general model with arbitrary lifetime distributions. 

Finally, in Section 4 we give the extinction results in our various different setups. Let us 
remark that, when the particle lifetimes have finite mean and the spatial dimension is small, 
local extinction of the population can be proved without the occupation times analysis; in 
this case a simple estimation yields the result, see the proof of Theorem 1. In contrast, the 
occupation time analysis is needed to treat the case of long-living particle types. 

2 Family tree analysis 

Following [3] (p. 553-558) we introduce the following auxiliary process. Consider a Markov 
renewal process with values in K, where in type i the process spends time according to a 



non-lattice distribution Fj (whose distribution function we denote again by Fj), such that 
Fj(0) = 0, and then jumps to type j with probabihty rriij. We write fii = f^ xTi{dx) for 
the mean of the i^^ hfetime, which can be infinite. Let tj{t) be the time that the process 
spends at state j up to time t. Put rij{t, a) = Pj {ij{t) > a}, i,j G K, where P, stands for 
the probabihty when the process starts in type i. We aim at finding an upper bound for the 
probabihties fij. 

First we show the connection between the Markov renewal process and the multitype 
branching system. We introduce the genealogical tree T of an individual, which comprises 
information on the individual's offspring genealogy, such as family relationships, mutations, 
death and birth times of individuals. For t > 0, let 7* denote the genealogical tree restricted 
to the time interval [0, t]. Finally, T{ stands for the reduced tree obtained from % by deleting 
the ancestry lines of those particles, which die before t. We write Pj for the law of T, if the 
process started from an ancestor of type i with age 0. From the context it will be always 
clear when P, refers to the branching particle system, or to the Markov renewal process. 

For any given t > and ancestry line w & %, let tj (w) > be the total time up to t that 
w spends in type j G K. Introduce the variable 

fij{t) = mintj{w), (2.1) 

which is the minimal time spent in type j among those particles that are alive at time t, 

with the usual convention that min0 = oo. We also define the maximum spent time in type 

j up to time t: 

aj{t) = maxtj{w). 
weTt 

(Notice that, in this case, the population procreated by the ancestor is not necessarily alive 
at time t). Let 

Uij{t,a) = Pi{fij{t) < a} 

denote the probability that starting from i, there is a particle at time t, who spent less than 
a time in j. Note that for t < a < oo, 

z/jj(t, a) = Pi{fij(t) < a} = Pi{ the process is not extinct at t} — )■ 0, 
as t — > oo, and for arbitrary a < oo, 

z/jj(t, a) = Pi{fj.j(t) < a} < Pj { the process is not extinct at t} — )■ 0, 
as t — )■ oo. Then by a renewal argument we obtain, for a < t, that 



z/,,,(t,a) = / F,(ds)[l-/,(l-z/.,j(t-s,a-s))] (2.2) 

Jo 

Ui,j{t, a) = 1-Ti{t)+ [ Ti{ds) [1 - /,(1 - u.j{t - s, a))] . 

Jo 



Since 1 — /i(l — z) < J2k=i ''^i,kZk, z = (2:1, ... , zk), we can compare the solution of (2.2) 
with the solution of the linear version 

na K 

ai^i{t,a) = / rj(ds)^mi,fcafc,i(t - s,a- s) (2.3) 

-^0 k=l 

aij{t,a) = l-Ti{t)+ rj(ds) V'mi,fcafcj(t- s,a). 

Notice that renewal argument implies again that <yij{t) = Pi{tj{t) < a} is the solution of 
the equation system (2.3). Let a\J{t,a) = i'ij{t,a), and let a'^"'^ = (ttj-"- )i,j=i,...,ft:, where 



afi (t, a) = Ti{ds) ^ rni^kdli {t- s,a- s) 

^0 k=i 

a^;'\t,a) = l-r,(t)+ / r,(ds)Vm,,ag(t 



By induction it is clear that Uij(t,a) < alj{t,a) for all n. We show that the iteration 
converges to the solution aij{t,a), and thus Vij{t,a) < aij{t,a). Let us fix a t > 0, and 
introduce the notation 

\\x-y\\t = sup{\xij{s,u)-yij{s,u)\ : i,jeK;0<u<s<t}. 
Then we get that, for all n, 



|a("+i)_ctW||, <max 

igK 



K 

t 



/ Vm,,fc||a(")-a("-^)||,r,(ds)=maxr,(t)||a(")-a("-^)||„ 



where we used that the mean matrix M satisfies Ylk=i "^«.fc ~ -'■ ^'^^ ^^^ ^- "^^^ ^^^^ estimation 
implies convergence to the solution of (2.3); the proof of uniqueness of solutions of (2.3) 
follows in the same way. Therefore we showed that z/jj(t,a) < aij{t,a). Notice that we 
only have shown that any solution of the equation system (2.2) is dominated by the unique 
solution of (2.3), which does not imply that (2.2) has a unique solution. We have proved: 

Lemma 1 For every a G (0,t) we have that 

Pi {3w e Tt : tj{w) <a} <Pi {ij{t) < a} , 

where the left side is for the branching process, while the right is for the Markov renewal 
process. 

Exactly the same way as in [3] Lemma 10, we can show a similar bound. 

Lemma 2 For every a G (0,t) we have that 

Pi {3weTt: tj{w) >a}<Pi {tj{t) > a} , 

where the left side is for the branching process, while the right is for the Markov renewal 
process. 



3 The Markov renewal process 

In this section we are going to analyze the auxihary Markov renewal process. First consider 
the discrete Markov chain Xi,X2, . . . with transition matrix M, and let p* = {pi, . . . ,Px) 
denote its stationary distribution. We have the following large deviation theorem for Markov 
chains ([5], Lemma 2.13.): For all 5 > there exist positive constants C,c, such that 



>6\ < Ce-'^", i G K, (3.4) 



n ■' 

where P, stands for the probability measure, when the chain starts from position i, and tj{n) 
is the number of visits to state j among the first n steps: 

t,(n) = #{/:X,=j,/ = l,2,...,n}. 

Here and below, several different constants arise in the calculations whose precise values 
are not relevant for our purposes. Therefore, for the reader's convenience we chose not to 
enumerate these constants. Hence the value of a constant may vary from line to line. In 
some proofs we use enumerated constants like ki,k2, ■ ■ ., whose values are fixed only in the 
corresponding proof. Finally, we use some global constants Ci,C2, . . ., whose values are the 
same in the whole paper. 

Let rit denote the number of renewals up to time t. With these notations we may write 



where Ci 1^2 i ■ • • ^"^^ ii^ random variables with common distribution function Fj, and 

^,U\ ^ \ ^ ~ ^nt, if Xrn = j, 

■' [^ 0, otherwise; 

that is Tij{t) is non-zero only for one term, and stands for the spent lifetime. Here Z„ is the 
sum of the lifetimes up to the n^^ renewal, and therefore Z„ is the sum of n independent, 
but not identically distributed random variables. 

3.1 A long living particle type 

Let 7 G (0, 1]. Assume that 

1 — Fi(a;) ~ x^"' , as a; — 7> 00 and (3.6) 

l-Tj{x)<Ax-^^, j = 2,3,...,i^, 

where A > and rjj > 1, j = 2,3, . . . ,K. Put rj = min{?7j : j = 2, 3, . . . , K}. We will show 
that with high probability the process spends ct times in type 1. 

Lemma 3 There exists Ci > such that for every ? G K and t > 1 

~^'^^^ < c,\ < Ct'-\ 



t 
for some C > 0. 



Proof. For simplicity we omit the lower index i. Recall that Ut stands for the number of 
renewals up to time t. For any A;2 > we may write 

The first term is easy to estimate. Due to (3.4), with probability > 1 — Ce^'^* we have 
ti{'nt)/nt > Pi/2, and so on this set 

^i(^) ^ ^h{nt)ti{nt) nt ^ '^tiK) k2P*i 
t ~ ti{nt) Tit t ~ ti{nt) 2 

Truncation and Cramer's large deviation theorem shows that for any d G (0, oo) there exist 
C, c > 0, such that for n Efi, 

" o(l) 

n 

Applying this with n ~ ct, the estimation above shows that the first term < C e"*^* (for some 
other pair of constants C, c) for any choice of Ci, A;2. 

Now let us investigate the second term. Clearly ti(t) < cit implies that tj{t) > k^t for 

some J > 2, with k^ = {1 - Ci)/K. If rit < k^t then t,(t) < 5|^J„^)+i < Sg.j+i by (3.5), 
therefore the probability in question is less then 



which proves our lemma. In the last step we used Theorem 2 of Nagaev [10], which says 
that for any c > and x > en, 

P {S^/'> -ni2j>x} < 2nx-"^A (3.7) 

for n large enough. In particular, for any 5 > 0, 



— ^ > 4 < cn^-'^^A. 



P 

n 



Notice that by using Lemma 6 below we obtain a stronger result. Namely, for any e > 

p|M^<Ci|<r+^-'?, 

for t large enough. Combining this with Lemma 1 we obtain 
Lemma 4 For any i G K there is a C\ and C > sueh that 

^<e,\< Ct'~^ 



3.2 Occupation times for j > 2 

To analyze the occupation times tj{t) for j = 2, . . . , K, we need a precise asymptotic for the 
number of renewals rit. 

We start by describing the asymptotic behavior of S'„ = ^i + ■ ■ ■ + ^„, the sum of n 
independent random variables with distribution function F, for which 1 — r(a;) ~ x~^, 
7 G (0,1]. In the following, limits of sequences are meant as n — )■ oo. We use the same 
convention for the continuous parameter t. 

Lemma 5 Assume that (i„ — )> oo if j < 1, and \ogn/dn — ;■ when 7=1. We have 

Moreover, for 7 < 1 there exist constants c^ such that for any sequences c^ for which c„ — )■ 
and n? ^^c„ -> 00, the following holds 

1 

1-7 
1/7' 



V{Sn< CnU^'"'] < 2exp 






// there exists a constant L > such that sup„n''' ^c„ < L, 7 G (0, 1], then for some c > 

P{5n<cX/^}<e-'=^ 

for n large enough. 

Note that in the case 7 = 1 we can choose c„ = c > arbitrary large. We will use this 
remark in the proof of Lemma 6. 

Proof. Let {t„} be a sequence of positive numbers such that n[l — r(t„)] — )■ if 7 < 1 
(i.e. tn = n^^^dn for some dn — )■ c>o), and nlogt„/t„ — )> when 7 = 1 (that is, t„ = ndn, 
where \ogn/dn — )■ 0). Using a theorem of Cline and Hsing ([1], Theorem 3.3) we get that 



lim sup 



'' s>t„ 



P{5n > S} ^ 



n[l-V{s)] 



0. 



This follows immediately from [1] if 7 < 1, while for 7 = 1 one has to check that the sequence 
which has to converge to 0, is 

n /■*" nlogtn 
— I xdi \x) ~ — . 



'^n J I 



Writing s = t„ in the form t„ = n}''^dn we obtain the first statement. 
Now we turn to the upper estimates for Snjn^^^ ■ 
We use a truncation method. For a > let denote 



' a, otherwise. 



8 



the truncated variable at a. For the first moment of this variable, as a — )■ oo we have 

^a = E^^^^ = rxdT{x) + a[l - T{a)] ~ ( p^'^"' f ' ^ ^ J' 
Jo I log«> for 7 = 1. 

For the proof we need Bernstein's inequality (cf. p. 855 of Shorack and Wellner [11]). Let 
Xi,X2, ... be iid random variables with EXi = 0, and let /t > and f > be constants 
such that E|X™| < vK,"^~'^m\/2. Then for the partial sum S'„ = Xi + ■ ■ ■ + X„ the following 
holds: 



P{|S„|>0<2exp{-^-iL_}. (3.8) 

Easy computations show that in our case (that is if 1 — r(a;) ~ x"'^) for m >2 

E[^('^)l™ ^^a'"-^ as a ^ oo. 

m — 7 

Since E|^*^") — fia]™' < E[,^*^")]™ + /i™, this shows that in Bernstein's inequality (3.8) we can 
choose V = 2a^~'^ and k = a. Obviously P {S'„ < c„n^/'''} < P ^ Sn < CnU^^'' k and so we 
may write 

P {Sn < CnU'/^} < P {Si^^ - rifla < C^U^'^ - n^^} < P {|5(") - nila\ > rifla - CnU'/^} , 

where in the last inequality we assumed that nfia > CnU^^"'- Applying Bernstein's inequality 
with V = 2 d^~'^ and /t = a we obtain 

P {Sn < Cnu'm < 2 exp I {nfi^-Cnn'/^Y \ 

Let nUa = 2n^^"'Cn, that is a ~ [2(1 — 7)]^^ n'^ Cn'"' =: a„. By our assumptions, a„ tends to 
00. Then the numerator in the exponential of (3.9) is n^^'^'c^, while the denominator is 

4([2(1 - -i)]-^ nht^f-^n + 2([2(1 - 7)]^ n^c^)c„n^/^ = c^n'^/^c^ , 
with c^ = (10 — 87) [2 — 27]!^. In this way we obtain finally that 

V{Sn<Cnn^'^] <2exp 






which is the desired bound. The last assertion in the lemma follows easily from Cramer's 
large deviation theorem, together with the truncation method. ■ 

Next we investigate the asymptotic behavior of the number of renewals rit in our Markov 
renewal process, where the lifetime distributions Fi, . . . , Tk are as in (3.6). We show that nt 
asymptotically behaves like the number of renewals fit in a standard renewal process, where 
the tail of the lifetime distribution is ~ x~"' . 

9 



Lemma 6 Let c(-) be a function such that 

lim c(t) = z/7 < 1, and lim c(t) logt = 2/7 = 1. (3.10) 

r/ien for all i E K, for t > /ar^'e enough, 

P.{^<c(t)}<2c(t), 0<7<1. 

7/7 < 1 t/ien /or any a G (0, 1 — 7), /or all i eK. and allt>0 large enough, 

If'~f= 1, ^/ien /or any a > and i G K i(/iere is a constant c > snc/i that for allt > large 
enough, 

Pj^>al<e-* 



' • t 



Proof. We drop the lower index i. In order to get an upper bound for rit let us define 

K 

f{x) = l[T,{x), 



j=i 



which is the distribution function of the lifetime ^= iaa,x{^^^\ . . . , ^^^^}, where ^^^\ . . . , ^^^"^ 
are independent and distributed as ri,...,rx respectively. Consider a standard renewal 
process n^, S'„ with this lifetime distribution. Recall that Z„, is the sum of the lifetimes up 
to the n^^ renewal. Clearly 

P {nt < t^c(t)} = P {Zirfc(t)\ > t} < P [Sirfc(t)\ > t} . 

According to Lemma 15 below, 1 — r(a;) ~ x~'^. Therefore using Lemma 5 we can write 

P{Vc«)J >t}<2c(t), 

where we have assumed that for 7 < 1 the convergence t/[t"'c{t)Y^'^ — )■ 00 holds, which is 
equivalent to c(t) — )■ 0, and that c(t) logt — > when 7 = 1. This proves the first statement. 

To obtain the lower bound we use the simple estimation Zn > 5*^ L-., that is, we simply 
drop the lifetimes with finite mean. First consider the case 7 < 1. Then we have 



{V+»j<t} C {<V-J)^^} 



^ S^[{pl-e)t'>+-i^^^ ^^T,+„J >Pl 2} \ [t^+^\ ^^ 2 



10 



hence using Lemma 5 and (3.4) we have 



^ P{Sm-s)t-'^^i^ '> + ''' 



o(l) 

— I / * \ ,., I „ 1 1 /_, / ~r c e 



L(Pi - £)i^+"J ^/^ ~ L(Pi - £)t^^''\ ^1^ 

< p<'_?LMzf)!:^!^<cr«/4+ce-^*"^" 

< 2exp|-ctA}+ce-"*"^". 

Taking into account that 7 + a > a/(l — 7), we obtain the statement. Finally, when 7 = 1 
we use exactly the same method. Using the event-decomposition as before, the last part of 
Lemma 5 and (3.4) we have 

P {^ > a} < P {5f,(,._,),j < t} + ce-^* < ce-^*, 

thus proving the last assertion of the lemma. ■ 

The preceding results allow us to obtain the following estimations for the probabilities of 
the smallness and largeness oitjit)/^. 

Lemma 7 Assume that (3.6) holds, and let c(-) he a function satisfying (3.10). Then for 
all i G K, any j = 2,3, . . . ,K, every e > and all t large enough, 

Pi {ij{t) < t^c{t)] <c{c{t) + f-^), < 7 < 1. 

// 7 < 1 then for any < a < 1 — 7, 

Pi {tj{t) > t^+"} < t^-"\ 

while for 7 = 1, 

Pi{tj{t) >at}< ct^^"^ 

for any a > 0. 

Proof. As before, for simplicity we omit the lower index i. Clearly, for any £ > 0, 
P {tj{t) < t"'c{t)] < P {tj{t) < t^c{t),nt > t^} + P K < t^} . 

Due to (3.4), on the set {rit > t^} we have, for any 5 > 0, that tj{nt)/nt G {p* — 5,p* + 5) 

with probability > 1 — ce~^^^ . Truncation method and Cramer's large deviation theorem 
show that 

11 



Since 

we have obtained that 

V{tj{t) <r'c{t),nt>f] <P{nt <ctM^)} + ce-'=*' < c (c(t) + e""*'), 

where in the last step we used Lemma 6. Taking into account that P {n^ < t^} < ct^~'^, 
which follows again from Lemma 6, the first inequality is proved. 

For the second part we use a similar technique. We first deal with the case 7 < 1. For 
any < 6 < 1 we may write 

P {tj{t) > t^+"} = P {tjit) > f'+\nt > t*} + P {tj{t) > t^^\nt < t^] . 

If n, < t\ using that f^^^ + ^,{t) < SJ^^l^^^, we get t,(t) < S^^^^^, < sgj^^. Therefore 
the tail probabilities of the right-hand term satisfy the inequality 

P {-^fiij+i > ^^^"} < 2^+''^ {[t^\ + 1) t-^^-^^+^M < t^~"^ (3.12) 

for all t large enough, where we used again Nagaev's result (3.7). Note that we only needed 
that 7 + a > 6. 

For the estimation of the other term we use again the decomposition (3.11). Since rit > t*, 
by (3.4) we have tj{nt)/nt G (p^/2, 2p*) with probability > 1 — ce~'^* . Due to Lemma 6, for 
any e' < a/2 

p{^>t^'}<2exp{-ctT^}. 

Since the orders of these terms are smaller than that of t^^'^'^ (given in the statement), we 
can work on {t^ < rit < T+^'j H {tjint)/nt E {p*/2,2p*)}. On this event, by (3.11) 



P - tj{nt) rit t-y - tj{nt) ^^ ' 



and so tj{t) > f^°- implies S^ ,.j^-^/tj{nt) > t" ^'/(2p*). Thus for t large enough 



>e 



t^ nt 

qij) 



where the last inequality follows again from (3.7). Choosing b such that a + 6 > 7 and 
e' < a + 6 — 7 we obtain the desired order t^"^"' . This, together with (3.12) gives the 
statement. 
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The proof in the case 7 = 1 follows a similar approach. For any 6 > 

P {tj{t) > at} < P {tjit) > at.nt < bt} +P{nt> bt} . 

The second summand in the right of the above inequality is exponentially small for any 
b > 0, and we have already shown in the proof of Lemma 3 that the first one is less than 
ct^^"^, provided that b is small enough. ■ 

Combining the last result with Lemma 1 and Lemma 2 respectively we obtain 

Lemma 8 Assume that (3.6) holds and that c(-) is a function satisfying (3.10). Then for 
all i G K, all j = 2,3, ... , K , each e > and all t large enough, 

Pi {3w G t; : tjiw) < rc(t)} = Pi {ijj{t) < rc(t)} < c (c(t) + 1'-^). 

If '-f < 1, then for any < a < 1 — 7 

P, {BweTf. t,iw) > t^+"} = P, {a,it) > t^+"} < t'-^\ 
while if J = 1, then for any a > 

Pi{3w G Tt : tj{w) >at} = Pi{cTj{t) > at} < ct^^'^. 

4 Extinction results 

In this final section, we apply the results on occupation times proved earlier in the paper to 
analyze extinction properties of our branching particle system. Let Nt denote the particle 
system at time t, i.e. A''^ is the point measure on R'^ x K determined by the positions and 
types of individuals alive at time t > 0. We write A^^* for the point measure representing the 
population of type-i particles at time t, that is Nl{A) = Nt{A x {i}) for any A C M^, hence 
Nt = N^ + ■ ■ ■ + N^ . As before the lower indices in P and E refer to the initial distribution. 
In particular, P^^i and 'E^^i refer to a population having an ancestor S(^x^i) of type i G K, 

initially at position x G M.'^. 

Let h : R'^ X K — )> [0, 00) be continuous function with compact support. We write 
(/i, h) = j hdfi for any measure fi on B(M'^ x K). Without danger of confusion we also write 

(x, y) = X]j=i ^iVi for the scalar product of vectors x = {xi, . . . , Xk) and y = {jji, . . . , y^)- 
Assume that the initial population A"o is a Poisson process with intensity measure A = 
Ai ^5{i} + ■ ■ ■ + \k^(^{k}, where i is rf-dimensional Lebesgue measure, and Aj, i G K, are 
non-negative constants. 

The Laplace transform of our branching process is, for any t > 0, given by 

E [e-<^"h)] = exp j - ^ A, f E,,, [l - e-^^"^^] dx 

= exp{-(A,l-E.,e-<^"h>)}. 
We put Ui{h,t,x) = E^,i (1 - e-^^*''^^). 
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To prove extinction of {Nf, t > 0} it suffices to show that the Laplace transform of A''^ 
converges to the Laplace transform of the empty population, and for this it is enough to 
verify that 

(A,f/.(h,t, ■)) -^ as t -> oo, 

which is the same as 

£/^(h.O:^/[E„(l-e-<«-.)]d.^O as < ^ oo fo. aU . e K. 

Let B G W^ he a ball, and assume that B x K. D supph. Then 

1 - e-^^*''^^ < I{Nt{B X K) > 0), 

which implies 

E.,, [1 - e-<^*''^>] < P,,, {NtiB X K) > 0} . 

Conversely, if hl^xK > 1? then 

1 _ Q-iNuh) y^^_ e-i)/(iVt(fi X K) > 0), 

and so 

E,,, [1 - e-^^^''^^] > (1 - e-i)P,,, {Nt{B x K) > 0} . 

In this way we get that 

Lemma 9 Extinction of {Nt, t > 0} occurs if, and only if for any bounded Borel set B C M.'^, 

/ Pi^^ {Nt{B X K) > 0} dx ^ for all ieK, as t ^ oo. 

Put a = min{aj : i = 1,2, . . . , K}. Recall the following result from [3]: 
Lemma 10 (Fleischmann & Vatutin). For each hounded B dW^ 

sup / E^_iiV((i? X K)dx — )■ as L t cx), 

t>l JR'i\C(t,L) 

where C{t, L) = {xeM."^ : \x\ < Lt^/'']. 

This means that extinction of {Nt, t > 0} occurs if, and only if for any bounded Borel 
set B (ZW^, and for L large enough 

Pi,x {Nt{B X K) > 0} dx ^ for alH G K as t ^ oo. (4.13) 

C{t,L) 
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4.1 Lifetimes with finite means 

When the hfetimes have finite mean and the dimension is small, it is not necessary to analyze 
the occupation times in order to prove local extinction. As we are going to show, in this case 
a simple estimation and the asymptotics of the extinction probabilities of critical multitype 
branching processes give the result. 

Let F^^' denote the probability generating function of the process starting from a single 
particle of type i: 

FW(t;si,...,Si^) = Ei s/' ■■■Sk , 0<Sj<l,jeK. (4.14) 

Put Q^'\t; Si,...,Sk) = 1- F^'^t; Si, . . . , Sk) and q^'\t; s) = Q^'\t; s, . . . , s). Clearly 

Pi,x {Nt{B X K) > 0} < Pj {the process is not extinct at time t} . 

Consider the discrete-time multitype Galton- Watson process {X„}, with the same offspring 
distributions as in the branching particle system. Let v and u respectively denote the left 
and right normed eigenvectors of the mean matrix M, which are determined by: 

vM = v, Mu = u, vu=l, lu = l. (4.15) 

Since by assumption M is stochastic, u = K^^l. Let f„ = (/^, . . . , f^) denote the generating 
function of the n^^ generation, that is /^(x) = Ej [x"^"] and put fi(x) = f(x). It is well- 
known that f„+i(a;) = f(f„(x)). Let us assume that 

X- {^r,l-f{l-ux)) ^x^+f^Lix) as a; ^ 0, (4.16) 

where /3 G (0, 1] and L is slowly varying at in the sense that lirax_^Q L{Xx) / L{x) = 1 for 
every A > 0. In this case, for the survival probabilities it is known that 

1 - f„(0) = (u + o{l))n-^/^Li{n) as n ^ oo, 

where Li is slowly varying at oo (see Theorem 1 in [12] or Theorem 1 in [13]). Moreover, 
assume that 

lim ;^\~ y ,;^j. = 0, ^ = 1,2,...,K. (4.17) 

n->oo (v, 1 - t„(0)) 

Then 

Q^''\t; 0) = Pj {the process is not extinct at t} ~ UiD'^t^^ Li{t) as t -)■ cx), (4.18) 

where D = '^^^lUiViHi; see Theorem 2 in [13]. 

Using the estimate above, we obtain the following theorem. 

Theorem 1 Assume that (4.16) and (4.17) hold. Then for d < a/f3 the process {Nf, t > 0} 
suffers local extinction. 
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Proof. Due to (4.18), for any e > 

Q^'\t;0)<ct-^, ieK. 
Plugging this into (4.13) we get 

Pi,^{A^t(5 X K) > 0}da; < ct^~^ . 

C{t,L) 

Since by assumption d < a/P, for some e > the exponent of t in the above inequahty is 
negative, which imphes that the integral in the left-hand side tends to 0. ■ 

Remark 1. When the generating functions /,, i G K, are of the form /j(s, . . . , s) = fi{s) = 
s + c (1 — sy~^'^^ where /3j G (0, 1], it is easy to verify that (4.16) holds with /3 = min{/3j : i G 
K}, and that (4.17) is fulfilled if for some e > 

lim n^+^+^[l - Ti{n)] =0, i = l,2,...,K. 

Remark 2. We remark that we do not need the precise asymptotic decay of the non- 
extinction probabilities given in (4.18); it suffices to know an asymptotic order of decay. In 
order to get this, instead of assuming in (4.16) that L(-) is slowly varying at 0, it is enough 
to suppose that L is an S-0 varying function, meaning that there exists an A > such that 
limsup^_^g L{Xx) / L{x) < A for any A > 0. S-0 varying functions were introduced by Drasin 
and Seneta [2]. The definition immediately implies that liminf^^o-^('^2:)/L(a;) > A~^ for all 
A > 0. It was shown in [2] that every S-0 varying function admits a representation as the 
product of a slowly varying function and a bounded (away from and oo) function. A careful 
analysis of the proof of Theorem 1 in [12] shows that, under the S-0 varying assumption on 
L, we have that for any e > and for all n large enough, 

|l-f„,(0)| <n-V. 



Since this estimate is precisely what we use in the proofs of our extinction theorems, all 
these results (including the infinite mean case) remain true in this more general setup. If 
(4.17) holds (which in particular implies that the lifetimes have finite mean), we obtain that 
for any e > and for all t large enough. 



Q'^'{t] 0) = Pj {the process is not extinct at t} < t 



l-E 

/3 



4.2 A lifetime with infinite mean — Case A 

From now on we assume that there is exactly one lifetime distribution with infinite mean; 
more precisely we assume (3.6). Moreover, in this subsection we additionally assume that 
a = min{aj : i G K} = ai, that is, the long-living particle type is the most mobile as well. 

In the following, P^ ■ denotes the distribution of the population starting with a single 
individual (5(^^j) of age 9 > 0. 
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Lemma 11 For all {x, i) eM.'^ x K, all bounded Borel B CM.'^ and all t > 0, 

Pi, {Nt{B X K) > 0} < C2 (t-'^/" + t'-") , 

where the constant C2 is independent of 6,x and i. 

Proof. Put A = {/Ui(t) < Cit}, i.e. A is the event that there exists a branch w E T{ 
such that ti{w) < C\ t, so the process spends less than ci t time in type 1 for some branch w. 
Clearly, due to Lemma 4, we may write 

P{Nt{B xK)>0} < P{A} + P{Nt{BxK)>0,A^} 

< ct^'"^ + E [Nt{B X K)Iac] . 

Conditioning on the reduced tree and noting that A is T/" measurable we have 

K Ni 



E^ [iV,(5 X K)/^c] = E^5^ Y. lArJ{W]{t)eB) 

j=i 1=1 






K Ni{Wi) 




)^ Y^ lA^KW^it) e B) t: 


i=i 1=1 


J 



E 



i=i 1=1 



where, given 7^' 



t ' 



riu\ 2. 



W{t)^W{t,,a,) + --- + W{tK,aK). 



(4.19) 



Here tj is the time that a branch of the reduced tree spent in type j, ti + ■ ■ ■ + t^ = t-, 
and {W{t^aj), t > 0} are independent symmetric aj-stable motions starting from 0, j = 
1, . . . ,K. Since on the complement of A any branch spent at least cit time in type 1, we 
have 

< ct-''/"^ = ct-''/'' 



(where pt-ti {x, dy) stands for pl^^ * ■ ■ ■*Pt^'^ {x, dy)), and we may continue writing the long 
equality as 



K 



< 



cr'^/'^Y^N}^ i^'^) <ct 



'd/a 



i=i 



Summarizing we obtain 



Pli{Nt{BxK) >0} <cr^/" + cti~''. 



Besides Lemma 11, our other key tool is an analogue of Lemma 3 in [15]. Recall the 
notations after (4.14). The proof is an easy multidimensional extension of the proof in [15]. 
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Lemma 12 If i] — I > d/a, then for any x E W^, t > 0, i G K and -u G (0, t — C2 ), 

P,,, {Nt{B X K) > 0} < g« {u; 1 - C2 (t - m)"'^/") , 
where the constant C2 is given in Lemma 11. 

Proof. Let \Nr\ = {N^{M.'^), ■■., iVf (M'^)), r > 0. For any u G (0, t), 

P,^,{Nt{BxK)>0} = Yl P.A\Nu\=k,Nt{BxK)>0} (4.20) 

k7^(0,...,0) 

= P.AINuIt^O}- Y1 PxA\Nu\=KNt{BxK) = 0}, 

k5^(0,...,0) 

where 
P,,a|iV„| = k,iVi(5 X K) = 0} = E [P {Nt{B X K) = 0||iV„| = k,&^,Y^} /(|iV„| = k)] . 

Here 6k is the vector of ages, and Y^ the vector of positions of individuals ahve at time u. 
Using independence and Lemma 11, the conditional probability inside the above expectation 
gives 



P{iV,(5xK) = O||iV„|=k,0k,lk} 
i=l 1=1 

i=i 1=1 

>{i-c,{t-ur'/i''\ 

where in the last inequality we used that 1 — C2 (t — ■u)^'^/" > 0. Therefore we obtain 

P,,,{|iV„| = k,iVi(fixK) = 0}> (l-C2(t-M)-'^/")"''p,{|iV«|=k}. 
Substituting this estimate back into (4.20), we finally get 

P.,,{NtiBxK)>0} < Yl [l-(l-C2(t-M)-"/")"'']pa|iV«l=k} 

keN'^vo 

= Q(^) {u; 1 - C2 (t - ix)-^/", . . . , 1 - C2 (t - u)-"/'') 
= g«(w;l-C2(t-M)"'^/"). 



Let us define the set 

A = {s G [0, If : f (s) > s}, (4.21) 



where an inequality of the form (xi, . . . ,xk) > {Vi, ■ ■ ■ .Vk) nieans here that x, > yi for 
1 = 1,2,. ..,K. 

We remark that, since 1 — f (1 — ux) < Mux = ux, we have 1 — ux G A for all x with 
< ux < 1. In our case u = K^^l, and this implies that the diagonal {(s, . . . , s) : s G [0, 1]} 
is contained in A. 

For given matrix families A{t) = {aij(t))ij and B{t) = {bij{t))ij, t > 0, let us define the 
matrix convolution C = A* B hj 



Cij 






The convolution of a matrix and a vector is defined analogously. Put Mp(t) = {mijTi{t))ij 
and recursively define 

M^+\t) = Ml{t)*M^^{t), n = l,2,.... 

Put also Mp(t) = (5jjr^(t))jj, where r°(t) is the distribution function of a constant random 
variable. Notice that Mp(t) constitutes the unit element in matrix convolution. The follow- 
ing multidimensional comparison lemma is borrowed from [13], which is a generalisation of 
Goldstein's comparison lemma [5]. 

Lemma 13 For any t > 0, any natural n and for all s G A, 

l-f„(s)-M;^*[(l-s)®r](t) < 1-F(t;s) 

n— 1 

< i-f„(s) + ^M^*[(i-s)®[i-r]](t). 

3=0 

Here x O y := (xi|/i, X2I/2, • • • , XKyx) if x = (xi, . . . , Xx) and y = {yi, . . . ^yx). 

We are going to use below the upper bound given in Lemma 13. The following lemma is 
Lemma 5 in [14]. 

Lemma 14 Consider two critical multitype branching processes sharing the same branching 
mechanism, with corresponding lifetime distributions r(t) = (ri(t), . . . , rx(t)) andr*{t) = 
(rt(t), . . . ,r^(t)). Assume that T{t) > T*{t) for all t > 0. Then for all t > and s G A, 

F(t;s)<F*(t;s), 

where F and F* are, respectively, the vector generating functions of the number of particles 
at time t in the first and second process. 

The main result in this section is the following theorem. 

Theorem 2 Assume that (4.16) holds, the mean matrix M is stochastic, and the lifetimes 
satisfy 1 — Ti{t) ~ t"'^ for some constant 7 < 1, and 

l-^,(x)<Ax-''^ j = 2,3,...,ir, 

where rjj > 1, j = 2,3, . . . ,K . Put rj = minjr^j : j = 2,3, . . . , K}. If rj — 1 > d/a and 
d < ^, then the process suffers local extinction. 

19 



Proof. Define the distribution function 



K 



f(t) = nr.(t), 



i=l 



wliicli is tlie distribution function of ^ = m.ax{C,i, ■ ■ ■ ,C,k}, where the random variables 
C,i, i = 1, . . . ,K, are independent with distribution function Fj. Lemma 15 below shows 
that 1 — r(t) ~ t^'^. Consider a new branching process where the branching mechanism 
is unchanged, but the lifetimes of all types have distribution F, and let F{t; s) denote its 

generating function at time t. Clearly, the choice of F shows that Lemma 14 is applicable, 
and so for s G A, 

F(t;s)<F(t;s). (4.22) 

(Notice that A, as defined in (4.21), depends only on the branching mechanism of our pro- 
cess). Now we apply the comparison lemma for this new process. Since now all the lifetimes 
have the same distribution, 

M?(t) = M"r"(t), 

where *" stands for the n-fold convolution. Moreover, for s = s 1, 
Ml * [(1 - s)(l - f)](t) = (1 - s){f*^{t) - T*^^+'\t))Mn = (1 - s){T*'{t) - T*^'+^\t))l, 

where we used the simple fact that M^ is stochastic if M is stochastic. Thus, in the rightmost 
inequality of Lemma 13 we get a telescopic sum, and therefore we obtain 

1 - Fit; Is) = Q(t; Is) < 1 - f„(ls) + (1 - s)[l - f™(t)]l. 

According to (4.16), for the survival probabilities we have 

l-/«(ls)<l-/«(0)<cn-i 
Taking into account (4.22) we have, for s G (0, 1), 

gW(t; Is) = 1 - F»(t; Is) < cn'^ + (1 - s)P {5„ > t} . 
Choosing n = t'^/(^+'^) and using that, by Lemma 5, 

P{Sn>t} = P (Sn > n^} < 2n-" = 2rT^/(i+"), 
we get 

g«(t; 1 - s) = g«(t; (1 - S)l) < ct'T^^ +st-^'/^^^'\ 

Hence, choosing u = t/2 in Lemma 12 we obtain the inequality 



q 



^^{u; 1 - C2 (t - m)"'^/") < criTW + crh-^. (4.23) 
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Multiplying by f^/", the second term in the right of (4.23) goes to 0, while in the first one 
the exponent of t becomes 

d 7 

a~ (l+e)/3' 

and this is negative ii d < a'-f/P and e is small enough. ■ 

The simple lemma we used above is the following: 

Lemma 15 Let X, Y be independent non-negative random variables with corresponding dis- 
tribution functions F and G. Assume that 1 — F{x) ~ x''^ and EiY < oo. Then for the 
distribution of Z = max{X, Y} we have 

l-H{z) := P {Z > z} r^ z-\ 

as z ^ oo. 

Proof. Since EF < oo, we have y{l — G{y)) — )> 0. Hence, 

z'^[l - H{z)] = z^[l-P {max{X, Y} < z}] 

= z^[l- F{z)G{z)] =z^l- F{z) + F{z){l - Giz))] 
= z^[l-F{z)] + z^F{z)[l-G{z)]^l. 



4.3 A lifetime with infinite mean — Case B 

Now let us investigate the case when ai is not the minimal a = min{aj : i = 1,2,..., K}. 
Without loss of generality, let us assume that a = a2. 

Notice that Lemma 11 is true in this case with exponent —d/ai, and so the variation of 
Lemma 12 also remains true. We state it for the easier reference. 

Lemma 16 If i] — 1 > d/ai, then for any x E ]R'^,t > 0,i G K and u E (0,t — Cg ), 

P,,, {NtiB X K) > 0} < g« {u; 1 - C2 (t - iz)-'^/"^) , 
where the constant C2 is given in Lemma 11. 
Put 

V = max J — , - I . (4.24) 

[ai a ) 

Lemma 17 Assume that 'yrj > d/a + 1. // 7 < 1, then for any e > 0, for any i E 
{1, 2, . . . , K} and for any bounded Borel set B, 



lim / P^i{A/'i(5 X K) >0}dx = 0. 

For 7 = 1 (then necessarily v = \Iol), 

lim lim sup / P^,i {Xt[B x K) > 0} dx = 0. 
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Proof. Without loss of generality, we will assume that 5 is a ball with radius r centered 
at the origin. First, consider the case 7 < 1. Put C{t) = {|a;| < t'"'^'^} and let e' < ae. Recall 
the definition of crj{t) after (2.1) and put 

A = {a^it) < t^+^', asit) < tT+^', ...,aKit)< t^+^'}, 

namely A is the set where, for all ancestry lines, the spent time in type j up to t is less than 
r+^' for all j = 2,3,...,K. 

First we work on the set A'^. By Lemma 8, 

K 

' -iri 



i=2 
According to Lemma 10, 

sup / Ea;,i [I{A'')Nt{B X K)] dx ^ as L ^ 00, 

t>l J\x\>Lt^/°' 

hence, it suffices to integrate on the region t"^^ < \x\ < Lt^^°'. On the other hand 

P,,,{^',Ar,(BxK)>0}<P{A'=}, 

and so 

f P^,i {A', Nt{B X K) > 0} dx < ct'^/^t^-"^ -^ 

due to our assumption. 

From now on we work on A. Translation invariance of the motion shows that 

/ E,,,/(A)iVi(B X K) dx = / Eo,./(A)iV,((5 - x) X K) da;. 

jR<'\C{t) JR''\C(t) 

By conditioning on the reduced tree, we can write 

Eo,^I{A)Nt{{B -x)xK) = Eo,/(iVt ^ 0)/(A) ^ J^ Pq, {wp{t) eB- x\r;} , 

j=l m=l 

where Wp(t) = W(ti, ai) + h Witx, o:k) as in (4.19). Integrating we obtain 

/ Eo,iI{A)NtiiB -x)xK) dx 

K NiiR-i) 

= Eo,/(iv, ^ o)/(A) J2J2 ^^ / Po'^ {^rw e '^yW'} ■ 

^•=1 ^=1 jRd\C{t) Jb^x 
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Since |a; + |/| < r and |a;| > t'"'^'^, we have \y\ > t"^^ — r > t"^^ /2 for all t large enough. Using 
Fubini's theorem and that j. k^*-^^ ~- '^(^) independently of y, the double integral can be 
bounded from above by 



c(r)Po, |W^rWI> 



ff+e 



V 



On the event A we can write 



4-v+e 



K 



^oA\WJ^m>^ T:\ < Y.v\\W{t,,au)\> 



fc=l 

K 



l-v+e 



2K 



V 



< J2^{\W{h,ak)\> 



<5 i-l/afc 



ft 



k=l 
K 



2K 



V 



fc=i *- ^ 



4.5 +l/a 



where we used that t'"'^'^ > t t^^ for some small enough 5 > 0, and that, by the definition of 

A and e', the inequalities t^^^ > t^t- > t^t- hold, while in the last step the self-similarity 
of the stable process was used. The last upper bound above goes to as t — ?> oo, and 



supE; 
oo 



K 



Nh 



< oo 



due to criticality of the branching. This finishes the proof of the lemma under the assumption 
that 7 < 1. The proof for the case 7 = 1 is a straightforward adaptation of the previous 
one. ■ 

The value 0:17 can be considered as the effective mobility of the type-1 particles. At an 
intuitive level if 0:17 > a, then second particle type is more mobile, even considering the 
long-living effect of the first one, so that in this case the "dominant" mobility is associated 
to the second particle type. The next two theorems deal with the cases when the first type 
is the dominant and when the second one, respectively. 

Theorem 3 Assume that (4.16) holds and that 'yrj > d/a + 1. If a > ai7, i.e. the mobility 
of the first particle type is dominant, then the process suffers local extinction for d < q;i7//3. 

Proof. Writing u = t/2 in Lemma 12, and proceeding in the same way as we did to obtain 
(4.23) in the proof of Theorem 2, we get 

q^'\t/2; 1 - ct-"^/"') < ct-'^/^ir^^/(^+^) + ct-^/(^+")^. 

Since in this case v = 1/ai, from Lemma 17 we get extinction provided that 



d 



< 



1 



ai (1 +£)/?' 
which holds for e small enough if c? < ai7//3. 
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Theorem 4 Assume that 'yrj > d/a + 1. //q;i7 > a, i.e. the mobility of the second particle 
type is the dominant one, then the process suffers local extinction for d < d^, where 

d^ = (mh 1 - (4.25) 

a Oil 

Proof. From the comparison lemma (Lemma 13) we have 

QW(t; Is) < cn'^ + (1 - s)P {Sn > t} . 

We have to choose t = n '< for some e > 0, and then minimize the estimations in e. In this 
case 

Putting u = t/2 in Lemma 16 renders 



j^'\t; 1 - s) < ct (1+^-)" +st'^. 



q^'\t/2; 1 - C2 1"'^/"^) < ct~(^W + ct"'^/"^-* 

Therefore we have to maximize 

. / 7 d , ^1 
mm < — -, 

\(l + e)/3'ai 1 + e 

with respect to e. Since the term 7/((l + e)f3) is monotone decreasing, and the term d/ai + 
£7/(1 + e) is increasing in e, easy computations show that the optimal choice is 

, 7(1 + r^) ^^ 

d/ai + 7 



and the estimation is 






Combining this with Lemma 17, and taking into account that v = '~f/a, we get extinction if 

a 1 + /3 

Solving the inequality, gives that extinction holds for d < d^, with the anticipated dimension 

d+. m 

Remark Notice that if 7/0; — l/ai — )■ 0, that is, if the effective mobilities of types 1 and 2 
are approximately the same, then d-^- — ?■ ai'-f/P, which is the critical dimension in Theorem 
3. Moreover, for fixed a, ai and 7, the critical dimension d^ considered as a function of /3, 
is decreasing, which is consistent with the known results. 

Acknowledgement. We are grateful to the referee for the comments and remarks that 
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